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Abstract 

We argue that ah evidences point towards a finite non-vanishing zero 
momentum renormahsed lattice gluon propagator in the infinite volume 
limit. We argue that different simulations with different lattice setups 
end-up with fairly compatible results for the gluon propagator at zero mo- 
mentum, with different positive slopes as a function of the inverse volume. 
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1 Introduction 

The lattice gluon propagator at small or vanishing momentum in the Landau 
gauge has recently been frequently addressed as it is related to several studies 
in the small momentum regime using non-lattice methods. It is often advocated 
that the zero momentum gluon propagator should vanish, while we have P shown 
a Slavnov- Taylor based argument in favor of a divergence when the momentum 
goes to zero. Notwithstanding these extraneous arguments we observe that the 
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genuine lattice data point towards a finite non-vanishing gluon propagator at zero 
momentum in the infinite volume limit. Our second claim is that, once a well 
defined renormalisation procedure has been defined, the different available results 
are close enough, despite several systematic effects, to suggest an agreement. 

Our aim in this note is simply to gather the arguments in favor of this claim, 
without discussing the relationship with any non-lattice claim. We will not 
present any new result but only quote published results and add a reanalysis 
of our old data. We concentrate on SU (3) pure Yang-Mills theory in the Landau 
gauge. 

There are two approaches to the problem. 

• One is to simply compute the gluon propagator at zero momentum and 
perform a well defined renormalisation. It is well known that the result is 
a finite non-vanishing value. But it might happen that the vanishing only 
happens in the infinite volume limit. Therefore an extrapolation to infinite 
volume is needed. 

• The second approach uses a set of small non-vanishing momenta and tries a 
fit of the propagator in terms of a power law or equivalently of the 
dressing function in terms of (p^)"'^. The fit gives some range of value for 
ac- The value ac = 1 - which corresponds to a no n- vanishing of the gluon 
propagator at vanishing momentum - has obviously zero measure and it is 
thus impossible to be assertive with this second method. It is nevertheless 
important to check that = 1 is compatible with the result and to check 
that the gluon propagator at vanishing momentum is in continuity with the 
result at small non- vanishing momenta. 

2 Definitions and notations 

In Landau gauge the gluon propagator writes 

G,Ap) = {5,u-^-f^)G^\p')^ (1) 

which implies 

G^'\V') = for PT^O 

G^'\0) = ^E^M.(0). (2) 

The factor 1/4 for zero momentum is due to an additional degree of freedom 
(clearly the orthogonality to the momentum in Landau gauge does not provide 
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any constraint for = 0) related to the fact that the Landau gauge fixing 
algorithm keeps unconstrained the global gauge transformations ^. 

In order to be able to compare the results from different gauge actions and 
different lattice spacings one needs to renormalise the gluon propagator. The 
standard method on the lattice is the Momentum substraction scheme (MOM) 
which amounts to define the renormalised propagator G^^^ from the bare one 
according to 

G(2)(p,a) = Z3(/i,a)Gg)(p,/i) (3) 
where the renormalisation condition is 

G^^\fi,fi) = whence Z3{fi,a) = /i^ G'-'^\fi,a) (4) 

A'' 

where a is the lattice spacing, i.e. the ultraviolet cut-off. 

This renormalisation can thus be done non perturbatively from lattice data 
provided that fi is in the available range for the given lattice spacing a. If this 
is not the case it is necessary to match the Za's with different lattice spacings. 
To illustrate this let us give an example: if we take fi = A GeV it is not possible 
to compute directly Z3 for the Wilson gauge action with jS = 6.0 ( = 1.97 
GeV) ^. We will thus use the results at /5 = 6.4 (a"^ = 3.58 GeV). We then 
need to compute Z^^fi, 6.0) /Z^i^fi, 6.4). This ratio is independent of fi at leading 
order. It can thus be computed non perturbatively for momenta in which both 
lattice spacings provide data. An analytic approach is to rely on the one loop 
perturbative formula 

Z3(/i,a') _ fP{a')Y'/'' 



Z3{ij,,a) \l3{a] 

which is valid for small enough lattice spacings (in the perturbative regime). 



3 The gluon propagator at vanishing momen- 
tum 

The Adelaide group has performed a systematic study [2] of the gluon propagator 
in the infinite volume limit. They use the mean- field (tadpole) improved version 
of the tree-level, O(a^) Symanzik improved gauge action. They choose a MOM 
renormalisation at /i = 4 GeV i.e. Gg^(4GeV,/x = 4 GeV) = 1/(4 GeV)^. They 
fit the volume dependence of the zero momentum gluon propagator on several 
lattice spacings and lattice volumes up to a volume of 2000 fm^, with always a 

•^Notice that this theoretically justified 3/4 factor is numerically confirmed as it ensures the 
continuity of the gluon propagator at — > which will be discussed later on. 
*It is advisable to keep p < {tt/2) a~^. 
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spatial cubic lattice and a length in the time direction twice the spatial length. 
Their fitting formula is 

Gf{0, ^ = 4 GeV) = GgL(0, = 4 GeV) + ^, (6) 

and gives 

Gg)^(0,/i = 4GeV) = 7.95±0.13GeV^^ c = 245 ± 22 fm^ GeV^l (7) 

This result clearly indicates a finite non vanishing G^^\o,fi). It is strange that 
nobody objects to this published result but that, nevertheless, one repeatedly 
reads that the zero momentum gluon propagator vanishes. 



p 


V in units of a 


bare propagator G*^^^ (p, a) 


1/L in GeV 


5.7 


16^ 


16.81 ±0.13 


0.0672 


5.7 


24^ 


15.06 ±0.29 


0.0448 


5.8 


16^ 


19.12 ±0.16 


0.0841 


5.9 


24^ 


18.12 ±0.30 


0.0685 


6.0 


32^ 


17.70 ±0.59 


0.0615 


6.0 


24^ 


19.67 ±0.35 


0.0821 



Table 1: Physical lattice sizes and raw data for the gluon propagator at zero 
momentum G^^\p, a) from our old data. 

This is why, waiting for a systematic and extensive reanalysis, we have simply 
digged out our old results for the gluon propagator which have been obtained 
from simulations with the Wilson pure gauge action on hypercubic lattices jSl El- 
Table^ lists the normalized raw data of the gluon propagator at zero momentum 
for our largest physical volumes (some of these data have never been published). 
No rescaling, perturbative (Eq. E)) or non-perturbative, has yet been applied 
to these data. Our volumes are not very large as this was not the aim of our 
simulations, and we do not claim our study to be an improvement over ref. |2] 
but simply an independent check. Using the same renormalisation as ref. [2] we 
find 

Ggj^(0,/i = 4GeV) = 9.1±0.2±0.2GeV"^ c = 140±30±40 fm^ GeV"^ (8) 

where the first error is statistical and the second is a systematic one estimated 
from different choices of the fitting points. 

More recently ref. [H] provides additional information on the same issue. In 
their table 2 the authors report fits of the zero momentum gluon propagator as 
a function of l/V, using only data at /3 = 6.0 with Wilson gauge action obtained 
on very anisotropic lattices ^. The results are given in lattice units and concern 
bare propagators. 

^The time length is typically 16 times the spatial one 
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reference 


Gg^(0,/i = 4GeV) in GeV-^ 


c in GeV ^ fm"^ 


max vol in fm^ 




7.95 ±0.13 


245 ± 22 


2000 


table ^ 


9.1 ±0.3 


140 ± 50 


90 


i 


10.9 - 11.3 


47- 65 


110 



Table 2: Summary of the infinite volume zero momentum propagator and its slope in terms 
of 1/V for three different simulations. The largest volume used in the fit is also indicated. The 
statistical error is not quoted in ref [S]. 



We shall assume, as has been done up to now, that the volume dependence 
is polynomial in 1 /V for large volumes. The very asymmetric shape is meant to 
provide very low values of the momentum ; it is interesting to check whether the 
zero momentum propagator depends on the geometry. We therefore convert the 
authors' fit of the 0-momentun propagator to physical units and perform a MOM 
renormalisation at 4 GeV for which we use a~^(/5 = 6.0) = 1.97 GeV and, from 
our non-perturbative fits : 

Zs{4: GeV, p = 6.0)) = 1.648, (9) 

and we get 

^2)^(0, /i = 4 GeV) = 11.3 GeV-2 and 10.9 GeV~^ 
c = 47fm^GeV-2 and 65fm^GeV-^ (10) 

where the two results correspond to a linear/quadratic fit in 1/V ^. We do not 
know the statistical errors. 

The results of these three collaborations are summarised in table. |21 Con- 
cernmg ^2)^(0, /i = 4 GeV) the three results are in the same ballpark and it may 
be conjectured that the systematic errors are not all taken into account: 0{a) 
effects, effect of the shape, insufficiently large volumes (for the second and third 
lines), uncertainty in the estimate of the lattice spacing in physical units, etc. 
Altogether it seems that, not only there is a clear indication in favor of a finite 
non vanishing zero momentum gluon propagator, but that different simulations 
agree on the value. Of course a more extensive study is necessary. 

Concerning the slope c the numbers clearly differ, they only agree in order of 
magnitude and are all positive. We expect that the slope is much more sensitive 
to systematic effects such as the shape. 

We turn now to the second approach, namely a fit of the p"^ dependence of 
the propagator at small momenta. We first claim that the gluon propagator is 
continuous and smooth at p = 0. This has been observed in several references 
(see for instance figure 17 in |2]) . This can also be seen in figure 2 in 

^We are aware that the authors of [Hj also use non polynomial fits in 1/V which can lead to 
vanishing or infinite zero momentum propagators. But they have themselves noticed that this 
destroys the smoothness. 
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The latter paper also compares the gluon propagator with periodic or twisted 
boundary conditions and concludes that the twisted propagtor is smaller than the 
periodic one but that the difference vanishes in the large volume limit. Let us now 
comment on the fit as a power law {G^'^\p) oc (p^)"'^^^) which necessarily discards 
the zero- momentum. In section 3.1 of ref. P we have shown with similar fitting 
formulae that ac is compatible with 1 on the examples of SU{2) and SU{3). But 
we have experienced instabilities and we do not know of any convincing results 
obtained with this method This instability may be due to the fact that, if such 
a power law applies in the small momentum limit, it can only be isolated at very 
small momenta which have not yet been reached. 

4 Conclusions 

The renormalised gluon propagator at zero momentum converges, in the infinite 
volume limit, towards a non vanishing finite value [2j if one uses a volume depen- 
dence which is polynomial in 1/V independently of the boundary conditions [7j. 
We have shown that different studies with different gauge actions, different pa- 
rameters and different shapes of the volume agree rather well on the value of the 
renormalised gluon propagator while the slopes in 1/^ agree only in sign and 
order of magnitude. 

No discontinuity of the gluon propagator is seen when the momentum goes to 
zero. It results that the infrared exponent ac (often named 2k) must be equal to 
1. We have stressed the instability of the fits of the propagator without the point 
at p = assuming a power law dependence: different fitting functions, which are 
equivalent in the infrared limit, give incompatible results. If this is duly taken 
into account in the systematic errors the value ac = 1 lies within the error bars 
in agreement with the claim about a non vanishing finite gluon propagator at 
zero momentum. 
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